Depth dependent dynamics in the hydration shell of a protein 
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We study the dynamics of hydration water/protein association in folded proteins, using lysozyme 
and myoglobin as examples. Extensive molecular dynamics simulations are performed to identify 
underlying mechanisms of the dynamical transition that corresponds to the onset of amplified atomic 
fluctuations in proteins. The number of water molecules within a cutoff distance of each residue 
scales linearly with protein depth index and is not affected by the local dynamics of the backbone. 
Keeping track of the water molecules within the cutoff sphere, we observe an effective residence 
time, scaling inversely with depth index at physiological temperatures while the diffusive escape is 
highly reduced below the transition. A depth independent orientational memory loss is obtained for 
the average dipole vector of the water molecules within the sphere when the protein is functional. 
While below the transition temperature, the solvent is in a glassy state, acting as a solid crust 
around the protein, inhibiting any large scale conformational fluctuations. At the transition, most 
of the hydration shell unfreezes and water molecules collectively make the protein more flexible. 

PACS numbers: 



I. INTRODUCTION 

Water is essential for the stability and function of pro- 
teins. It modifies properties of proteins such as the dis- 
play of the dynamic transition temperature, T^. The 
dynamical transition corresponds to the glass transition 
observed in the hydration shell of a protein accompa- 
nied by a correlated rubber-glass transition of the pro- 
tein. The transition results in enhanced atomic fluctu- 
ations with respect to dehydrated proteins and conse- 
quentely increased flexibility of the backbone.— The dy- 
namical transition temperature depends on the time scale 
of the observation. Since the observed configurational 
fluctuations are delimited by the resolution of the mea- 
surement, the reported value of the dynamical transition 
temperature depends on that resolution.^ For example, 
neutron scattering data yields Td = 240 K at a time 
window of 50 ps while it decreases to Td = 220 K at a 
time scale of 2 ns.— >^ The transition is accompanied by an 
increase in the slope of the mean squared atomic fluctua- 
tions. The local relaxation times of Ca atoms is observed 
to change in accord with the increased flexibility of the 
backbone.- However, in the absence of solvent, the dy- 
namical transition is not observed^ and to recover the 
transition the charged residues of the protein have to be 
covered by solvent 

Proteins in turn affect the dynamics of water in the 
hydration shell. Experiments at room temperature show 
that the relaxation times of water molecules in the vicin- 
ity of proteins are larger than in the bulk^ and highly de- 
pendent on the location around the protein.— While cou- 
pling of water dynamics with the onset of the dynamical 
transition has been investigated^itiS, and it is widely ac- 
cepted that the solvent dictates the transition, the mech- 
anism by which water operates on proteins is still an open 
question.^"* Recently, single molecule rotational correla- 
tion time of water molecules around the protein has been 
determined from spin relaxation experiments at various 
temperatures. This is more localized than translational 



diffusion, and in turn is a more accurate measure of mo- 
bility in the hydration layer. Since it is crucial to have 
a better understanding of the dynamics of water and its 
effect on the protein both below and above T^, we study 
the translational and rotational motion of water in the 
vicinity of the protein in these separate regions, along 
with that of pure water and the protein. This study en- 
ables us to understand water dynamics on the time scales 
that are relevant to protein motions. 

The paper is organized as follows; In Sec. [Ill we 
describe the numerical methods and the details of the 
molecular dynamics (MD) simulations carried out for 
pure water and proteins in water. In Sec. IIII[ we flrst 
determine the glass transition temperature of pure water, 
described by the TIP3P model. Thereafter, we investi- 
gate the local dynamics of two commonly studied pro- 
teins, lysozyme and myoglobin, at room temperature as 
well as at 180 K. The latter is below the dynamical tran- 
sition temperature of common proteins, but above the 
glass transition of water. The hydration levels and water 
residence times around the residues are computed, and 
the predominant contributions to the observed behavior 
are shown to be accounted for by two simple models, de- 
veloped in this work. Afterwards, the Ca and water re- 
laxation times and their temperature dependence are in- 
vestigated to uncover the nature of the coupling between 
local protein dynamics and water. Finally, conclusions 
are drawn in Sec. IIVI 



II. NUMERICAL METHODS 

To study the dynamics of water in the hydration shell 
of a protein, we carried out MD simulations for two dif- 
ferent proteins, lysozyme of 129 residues and apomyo- 
globin of 152 residues (Protein Data Bank codes 61zm 
and ljp6, respectively). Observing common patterns of 
a hydration shell for these two proteins, one may de- 
duce universal properties of hydration water dynamics. 
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The proteins are solvated so that the water to protein 
mass ratio is h = 3.42 and h = 3.34 for lysozyme and 
myoglobin, respectively. These values are we above the 
minimum of h = 0.7 for which a fully hydrated protein is 
observed, e.g., for bovine pancreatic trypsin inhibitor.— 

The solvated simulation boxes are equilibrated for 2 ns 
at 300 K, and then a further 2 ns at the desired tempera- 
ture. All simulations are in the NPT ensemble at a pres- 
sure of 1 atm. Two sets of simulations are carried out, 
one at 180 K, below the dynamical transition tempera- 
ture, and the other at 300 K. The NAMD packageis is 
used with the CHARMM 27ii force field and the TIP3P 
water model. A time step of 2 fs is used and 50 ns long 
trajectories are produced. Coordinates are recorded ev- 
ery 0.4 ps. Further details on the simulation methods are 
as in Okan et al}^ 

We investigate the equilibrium characteristics and the 
dynamics of water molecules in the vicinity of each 
residue. The hydration shell around a given residue is 
defined to contain the water molecules within 6A of its 
Ca atom. Smaller values of the cutoff results in a too 
small number of solvent molecules, and hence poor statis- 
tics, while larger cutoff radius will result in taking into 
account more of the bulk dynamics. We note that the 
correlations in the solvent can extend up to a distance of 
lOA away from the protein surfaceJ^i^ 

To differentiate the relative importance of the bulk ver- 
sus the local protein environment, we display the char- 
acteristic properties of the residues as a function of the 
depth index, d?^ This index is defined as twice the ratio 
of the accessible volume of an atom Va to the exposed 
volume of an isolated atom Vq, 



Thus, the deeper the residue is in the protein, the smaller 
is the depth index. A surface residue has approximately 
half the accessible volume of an isolated atom, and its 
depth index is close to one. For an isolated atom in bulk 
water, the depth index is two. The exposed volume Va is 
calculated in a sphere of radius r. If the radius is chosen 
too small or too large, all values of depth indices converge 
to or 2. The optimal value of r corresponds to the case 
when only one residue has a vanishing depth index; for 
proteins in this work one finds r = 9A. In all the results 
presented as a function of depth index, we average over 
residues with similar depth indices, where the bin sizes 
are evaluated from dmax/lS. The resulting error bars are 
shown on the figures. 

To evaluate the temperature dependence of the dynam- 
ics we also perform MD simulations of TIP3P bulk water, 
and hydrated lysozyme for a series of temperatures. For 
the former, we span the temperature range of 80 - 350 K 
at 1 atm, for 10 ns on a system of 273 water molecules. 
For the latter, we perform the MD simulations in the 
range 160 K to 300 K, each of length 24 ns, keeping all 
other conditions the same as the protein-water systems 
described above. 



III. RESULTS AND DISCUSSION 

A. Hyration levels and water residence times 

One can measure the hydration levels of each residue, 
defined as the average number of water molecules within 
QA of its Ca atom present. The number of water 
molecules in the hydration shell of a residue is highly 
dependent on depth: The deeper inside the protein, the 
higher the coordination number of a residue is and the 
lower its accessible volume. Although the solvent density 
p depends on depth, i.e. the first layer around the pro- 
tein has a density about 15% larger than bulk water 
we assume p does not depend on depth index as a first 
approximation. The number of water molecules is 
then proportional to the accessible volume, Vq d/2, 

N^ = ^d. (2) 

Since the bulk density p is known, one may predict the 
number of water molecules as a function of the depth 
index using this relationship. We compute the average 




Depth Index 

FIG. 1: Number of water molecules as a function of the depth 
index at 180 K and 300 K for lysozyme and myoglobin. The 
straight lines are obtained from Eq. (2] 

number of water molecules both in the hydration shell 
of each residue and in the bulk. The bulk values are 
32.4 and 30.5 water molecules at 180 K and 300 K, re- 
spectively, corresponding to a density of 1070 kg/m'^ and 
1007 kg/m^. The results are depicted in Fig. [1] and are 
in agreement with the linear dependence predicted by 
Eq. [5] for both proteins. Hence, the temperature depen- 
dence of the distribution of water molecules in the protein 
is determined predominantly from the bulk value p{T). 
Local structural irregularities and protein dynamics have 
a relatively minor role. 
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FIG. 2; Normalized number of water molecules initially re- 
siding in the vicinity of a residue as a function of time. Sam- 
ples are for residues with different depth indices (CYS30: 
d = 0.065, ASN65: d = 0.435, THR47: d = 1.26) and in 
the bulk (d — 2). The dashed lines represent the stretched 
exponential fits, d increases from the top towards the bottom 
curves. 



FIG. 3: Residence times of water molecules as a function of 
the depth index for lysozyme and myoglobin from fits to Eq. 
O The stretched exponents fluctuate around 0.5 at all depths 
and temperatures for both proteins. The limiting value in the 
bulk is represented by the dashed lines. Predictions from Eq. 
[S] (solid line) and from Bezrukov el al— (dot-dashed line) are 
also shown. 



Although the distribution of water molecules, an equi- 
hbrium property, is not affected by temperature, the dy- 
namics of water around the protein depends strongly on 
T. We characterize the latter by their residence times, 
Tr^ which is the average time a molecule takes to escape 
from a given region. To measure the residence times of 
water molecules in the hydration shells of all residues, 
we partition the 50 ns long trajectories into fifty 1 ns 
long pieces, and record the initial number of hydration 
shell water molecules in each region. We then compute 
the decrease in the number of these water molecules as a 
function of time by checking if they remain in the hydra- 
tion shell. Since the trajectories are recorded in intervals 
of 0.4 ps, water molecules that cross the boundary and 
return to the monitored layer faster than that time scale 
are assumed to have remained. We depict in Fig. [5] this 
decrease for selected residues of lysozyme, a highly buried 
one CYS30 [d = 0.065), one at an intermediate depth 
ASN65 {d = 0.435), a surface residue THR47 {d = 1.26) 
and finally the escape process in the bulk {d — 2). Many 
different processes contribute to the curves in Fig. [21 
Most of the water molecules are close to the edge of the 
sphere we consider, and thus can escape faster on av- 
erage. On the other hand, the water molecules at the 
center of the sphere initially will remain longer. There 
will also be additional contributions from the fluctuating 
chain molecule. Consequently, one expects to have a su- 
perposition of exponential decays for N^it) which is best 
described by a stretched exponential function;^ 

N^it) = {N^) eicp[-it/Trf], (3) 



where is the effective residence time and /? is the 
stretched exponent. The fits are illustrated in Fig. [2]and 
are in good agreement with the MD results, validating 
the model. The stretched exponent has a value /3 w 0.5 
regardless of the location of the residue and temperature. 
We thus compute the residence times for all the residues 
of lysozyme and myoglobin, and plot them as a func- 
tion of depth index in Fig. [31 The residence times for 
equivalent depths can have important variations because 
of the residue specificities and their spatial arrangements. 
For example, residues can be hydrophobic or hydrophilic, 
and consequently decrease or increase the residence time 
of the surrounding water. Moreover, residues can be con- 
centrated in one region of the hydration shell facilitating 
the escape process. However, one can capture the general 
trend when considering a large number of residues. Dif- 
ferences of up to three orders of magnitude between the 
residence times of solvent molecules buried in the protein 
compared to those near the surface are observed. This is 
due to the fact that the water molecules inside have to 
escape traps formed by the residues. Thus, the higher the 
coordination number of a residue, the larger the number 
of traps, and the longer the escape time. Note that, for 
certain residues that display extremely slow decay pro- 
files at 180 K (exemplified by the uppermost curve in 
Fig. [2]), Eq. [31 provides residence times longer than 100 
ns. This is due to the error involved in extrapolation of 
the initial decay profiles to very long time scales. 

Using a model of diffusive motion in the presence of 
local obstacles, a relationship may be derived for the res- 
idence time as a function of depth. As a first order ap- 
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proximation, we assume the probability to escape from 
the neighborhood of a residue Pr decreases hnearly with 
the occupied volume as 



P, = Po 1 



Vo 



(4) 



where Pq is the probability of escaping in the bulk and 
Vr = ^ occupied volume around the residue. 

The occupied volume is the total volume minus the avail- 
able volume Vr ^ Vq - Va, using d = 2Va/Vo (Eq. [1]), 
one gets Pr = Pq d/2. Since the residence time scales in- 
versely with the probability of escapCj^- and the propor- 
tionality constant is assumed to be an intrinsic property 
of the escaping species irrespective of the local environ- 
ment, one gets, 



Tr{d) 



(5) 



where tb is the residence time in the bulk. In Fig. [31 
we display the findings from Eq. [S] along with the results 
from the highly elaborate analytical solution of Bezrukov 
et al?^ for the probability of diffusion displacement of a 
molecule using the fraction of possible trajectories in the 
presence of obstacles. The rather simple model of Eq. 
[S] captures the main features of this detailed analytical 
solution. We further find that both approaches predict 
the depth dependence at 300 K, while they fail to do so 
at 180 K. 

So far as the predictions of Eq. [S] reproduce the main 
features of water behavior around the protein at room 
temperature, they corroborate that specificity plays a 
lesser role in determining residence times at physiolog- 
ical temperaturesj^i^ Below T^, however, we find that 
the simple trapped diffusion idea does not capture the dy- 
namics of water. In fact, the calculated residence times 
are approximately one order of magnitude slower than 
predicted by the theory, indicating that the dynamics is 
not controlled by diffusion in the presence of obstacles, 
hence suggesting that hydration shell water is not liquid 
below the transition temperature. This confirms the ob- 
servation that translational diffusion is hindered below 
the dynamical transitioui^^ 



B. Depth dependent relaxation times in the 
hydration layer 

To further investigate the water/protein interactions, 
we quantify the dynamics of the protein, hydration layer, 
and bulk solvent by measuring relaxation times. For the 
Ca atom of the ith residue, we follow the decay of the 
correlation function, 



(AR,(0-AR,(0)) 
(AR|) 



(6) 



where ARi(<) = Ri(i) — (R^) is the fluctuation vector of 
the atom. We note that (R^) is computed for each 0.5 
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FIG. 4: Ca correlation function for residues with different 
depth indices. The circles are for the buried residue (CYS30, 
d — 0.065), the squares for the intermediate (ASN65, d — 
0.435) and the diamonds for the surface one (THR47, d = 
1.26). The straight lines represent the stretched exponential 
fits. 



ns trajectory piece and then a best fit superposition is 
performed to separate the internal motion of the protein 
from long time tumbling. Previously, this method per- 
mitted evaluating the dynamical transition temperature 
for different protein-solvent systems i^i^ 

Many different processes over a wide range of time 
scales contribute to the backbone dynamics as we re- 
cently discussed in detail Here we employ stretched ex- 
ponential functions on the correlation functions of Eq. [6] 
to fit the initial decay processes as in our previous work.— 
We depict in Fig. Uthe correlation function for the three 
residues of lysozyme also exemplified in Fig. [2] and the 
corresponding stretched exponential fits. We note that 
the long time tails of these relaxation functions belong to 
more collective dynamics with a nanosecond time scale, 
which we do not further elaborate upon since this was 
discussed in great detail recently.— Closely associated 
with backbone dynamics, we further compute the water 
relaxation times by using the local polarization vector 
in the vicinity of each Ca as depicted in Fig. [S] One can 
then compute the following correlation function. 



cnt) = 



{AP,{t)AP,{0)) 



(7) 



where APi(i) = |P^(t)|-(|Pj(t)|) is the fluctuation of the 
norm of the polarization vector of the ith residue. The 
polarization vector is the sum of the polarization vectors 
of the water molecules within the hydration shell of the 
residue considered, 



(8) 
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FIG. 5: Schematic representation of the hydration shell of . 
residue and the local polarization vector P. 



We emphasize that we use a difference definition utilizing 
the norms ([7]) , as opposed to the usual expression involv- 
ing the scalar product of the dipolar moments. While 
the dipolar moment in the bulk is equal to zero on av- 
erage, it is not the case inside the protein, where water 
molecules can have preferred orientations because of the 
proximity of protein atoms. Consequently we use the 
quantity APi{t) which eliminates this bias, in addition to 
the advantage of Eq. ([7]) providing a faster decay since 
it corresponds to a higher moment. Moreover, we note 
that the relaxation time obtained depends on the num- 
ber of water molecules considered for either definition of 
the relaxationi^ To recover the experimental value of the 
relaxation time, one should compute the polarization vec- 
tor for larger spheres. At room temperature, we find that 
the radius should be larger than the correlation length of 
the solvent which is ca. 12 A. Beyond this sphere size, 
the Debye relaxation time converges to the value of 7.3 ± 
0.7 ps, as also reported in literature previously for TIP3P 
water model-^ , and compares well with the experimental 
value of 8.2 ps. In contrast, for the sphere size chosen 
in this work, 6 A, the value is 5 ps. For the relaxation 
function of Eq. ([7]), the respective relaxation times are 
3.2 and 1.8 ps for spheres of radius 12 A and 6 A, respec- 
tively. Thus, our computed values of the relaxations are 
consistent within this work, and are useful for comparing 
the local dynamics of water clusters at various depths. 
However, they cannot be compared directly with Debye 
relaxation times. As for the backbone, we use stretched 
exponential fits to the initial decay region of the correla- 
tion function to extract the relaxation times. We depict 
in Fig. ini for the three typical residues and the bulk, the 
correlation functions along with the corresponding fits. 
We note that, here we study the dynamics of a group of 
water molecules, and not the dipolar relaxation of indi- 
vidual water molecules. The former permits us to study 
collectivity and hence evaluate the glassy behavior of wa- 
ter clusters around individual residues. 

We observe that for both the backbone and the wa- 
ter molecules, the short time relaxation is relatively well 
described by the stretched exponential functions. For 
longer times a slower, single exponential decay also ap- 
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FIG. 6: Local water polarization correlation function as a 
function of time for different depth indices. The straight lines 
represent the stretched exponential fits. Symbols are the same 
as in Fig. [l] Triangles are for bulk water. 
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FIG. 7: Relaxation times of Ca atoms and water molecules 
as a function of the depth index for lysozyme (circles) and 
myoglobin (squares). The straight lines indicate bulk water 
relaxation times. 



pears in the water correlation functions, attributed to the 
large scale reconfiguration of the water molecules inside 
the region. However, for the current purpose of studying 
water/backbone dynamics on the picosecond time scales, 
we focus our attention to the former relaxations of both 
the backbone and the hydration shell water. Having val- 
idated the model for the correlation functions, we com- 
pute the relaxation times for each residue of lysozyme 
and myoglobin. We depict them as a function of depth 
index in Fig. [7] and their corresponding stretched expo- 
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FIG. 8: Stretched exponents of Ca atoms and water molecules 
as a function of the depth index for lysozyme (circles) and 
myoglobin (squares). The straight lines indicate bulk water 
values. 



nents in Fig. |5] for the two temperatures. We observe 
that the relaxation times of Cq atoms (tq.) are higher at 
room temperature than at 180 K (the average values are 
13 ps and 6 ps, respectively). We have previously ex- 
plained this slowing down of the local relaxations using 
a model that relates the onset of the coupling between 
the local internal motions of the protein to the hydration 
layer with temperature. The results showed that such 
behavior is due to the interplay between the decreased 
stiffness and the modified effective friction coefficient.— 
Moreover, independent of temperature, increases with 
depth index, which may be explained by packing argu- 
ments: The more buried a residue is, the larger is the 
coordination number, and hence the local rigidity of the 
medium. Consequently, the magnitude of fluctuations is 
smaller, leading to faster relaxations. 

Conversely, relaxations of water molecules forming the 
hydration shell (r^,) are observed to speed up with in- 
creasing temperature. At 180 K, their values are highly 
depth index dependent, accompanied by slowly increas- 
ing stretched exponents with values close to /? — 0.2. The 
latter is indicative of a collection of trapped states for 
the rotational motion of water molecules according to a 
model of hierarchically ordered materialsi^"— Thus, the 
water layer shows a strong tendency to be affected by 
the local environment (i.e. packing) below T^. However, 
their dynamics is independent of depth index at 300 K 
beyond a critical depth index (d = 0.3), where there are 
at least five water molecules surrounding a given residue 
(see Fig. [T|), marginally enough to form a polarized "wa- 
ter layer". The relaxations have a characteristic time 
scale of approximately 5 ps at room temperature. The 
distribution of relaxation times is much narrower as in- 
dicated by the larger stretched exponents of /3 > 0.4 (see 



Okan et alJ^ for a detailed interpretation of stretched 
exponents). 

The fact that r^, are on the same order as Tq, at 300 K, 
along with their depth independence suggests that water 
and overall protein dynamics are coupled at physiological 
temperatures, supporting the idea that water forms a 
vicinal layer around the protein,^ accompanied by a large 
water reorganization energy,^ In contrast, r^, values are 
1 — 2 orders of magnitude slower than Tq at 180 K, along 
with a very strong dependence on the local environment. 

Moreover, at room temperature hydration shell water 
relaxes slower than bulk water, compatible with a more 
viscous layer, slowed down by the interactions with the 
residues. On the other hand, below the transition, due 
to the lack of larger displacement within the hydration 
shell, hydration water relaxes faster than the bulk. This 
suggests that within the time scale of the observations 
leading to the correlation function, the hydration shell 
water is in a glassy state. 

These findings corroborate the results on the residence 
times in Fig. [31 They lead to a better understanding of 
the relationship between water-protein interactions and 
the dynamical transition. We further investigate this by 
computing the relaxation times of the motions described 
by Eq. [7] for lysozyme over a wide temperature range. 



C. Gradual unfreezing of the hydration layer 

We have produced 24 ns trajectories for lysozyme at 13 
separate temperatures spanning 160 to 300 K and com- 
puted the relaxation times of the hydration shell and bulk 
water. We observed in Fig. [7] that below T^, hydration 
shell water relaxes faster than the bulk if the residue's 
location is deep enough. To observe the transition from 
faster to slower relaxation, or in other words from glassy 
to liquid, one can simply count the number of residues 
slower than the bulk at each temperature as depicted in 
Fig. ini One can then measure the critical depth index dc 
beyond which the hydration shell water becomes slower 
than the bulk, as shown in the inset to Fig. [SI We observe 
a marked transition from a state where all the hydration 
shell waters are faster than the bulk to the opposite. The 
observed behavior may be described by a gradual unfreez- 
ing of the hydration shell. At temperatures just above 
Tg of water (see the Appendix for a prediction of the 
Tg of the TIP3P water model in the range of 150 - 160 
K), only the outermost hydration shells are not frozen. 
For example at 180 K, residues of depth index approxi- 
mately larger than 0.5 have equivalent relaxation times 
to bulk water (also see Fig. [7]). We note that some wa- 
ter molecules (~10%) remain faster than the bulk even at 
temperature well-above Td] these are highly buried water 
molecules that remain coupled to the protein. The grad- 
ual unfreezing of the hydration shell enhances our un- 
derstanding of the dynamical transition occurring within 
a wide time and temperature range, determined by the 
time scale of solvent coupled structural relaxations. It 
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cally relevant motions. 

Bulk and hydration layer water have been shown to 
separately control motions in functional proteins; for ex- 
ample, the former allows ligand entry/exit, while the 
latter is responsible for migration of ligands within the 
proteini^ The hydration layer acts as a lubricant for 
the onset of the functional dynamics in the protein, as 
new channels between conformational substates emerge, 
e.g. through jumps enabled in the side-chain torsional 
angles. '^^ Our results demonstrate that hydration water 
not only performs localized motions, but also participates 
in the global dynamics through long-range diffusion. Be- 
low the transition temperature, on the other hand, the 
miscommunication between the substates maintains the 
trapped hydration water highly dependent on the dynam- 
ics of the protein. 



FIG. 9: Normalized number of residue for which the hydra- 
tion shell relaxation time is slower than that of bulk water. 
The dashed line represents a Boltzmann sigmoid fit yielding 
a midpoint value at 187 K. The inset represents the critical 
depth index dc as a function of temperature, the line is a guide 
for the eye. 



has already been suggested that the broadening in the 
glass transition of the protein hydration shell is due to a 
distribution of water clusters with different glass transi- 
tion temperatures."^^ The heterogeneity of the hydration 
shell quantified by the depth index and its relation to the 
hydration levels, residence times and relaxation times ex- 
pands on this view. 



IV. CONCLUSION 

We study the residence times and memory loss of water 
molecules around individual residues to interpret their 
role in the dynamics of folded proteins. The calcula- 
tions are repeated for two well-studied proteins, leading 
to identical results: Above the residence time of water 
molecules around a given residue may be interpreted via a 
simple depth dependent model of trapped diffusion. The 
memory loss of the polarization vector displays a simi- 
lar trend to that of the protein backbone in terms of the 
time scales involved. In contrast, the dynamics of wa- 
ter molecules are non-diffusive in the regime below Td- 
Consequently, below Td the hydration shell is solid-like, 
behaving as a crust around the protein. Upon heating, 
the hydration shell is unfrozen, and permits the transfer 
of entropy from the bulk to the protein by the entry / exit 
process of water molecules, hence acting as a plasticizcr 
and increasing the overall flexibility of the protein. The 
contrast between the water dynamics in the hydration 
shell of a protein at physiological and non-physiological 
temperatures therefore reveals the interactions and the 
free energetic requirements necessary to achieve biologi- 



V. APPENDIX 
A. Glass transition temperature of TIP3P water 

One of the crucial quantities regarding the dynamical 
transition in proteins is the water glass transition tem- 
perature. The glass transition temperature is defined as 
the temperature at which the shear viscosity is 10^^ Pa.s, 
corresponding to a relaxation time of approximately 100 
s. It is also marked by a crossover from ergodic to non- 
ergodic behavior. Experimentally, a step-like decrease in 
the susceptibilities (e.g. heat capacity) is observed. It 
was suggested that the dynamical transition of the pro- 
tein is dictated by the glass transition of the solvent i^"— 
Although pure water readily forms the ice-phase and is 
very difficult to vitrify, an experimental value of ~ 165 K 
has nevertheless been measured j '^^'^'^ To evaluate the in- 
fluence of glassy water on the protein, we must first assess 
the behavior of the water model used in this study. 

The specific heat may be computed from 



and the isothermal compressibility as 



ksTiV) 



(9) 



(10) 



where E is the total energy of the system and V the 
volume. The results along with experimental values are 
depicted in Fig. (TU] The specific heat and compressibility 
curves yield a glass transition temperature between 150 
K and 160 K. This value is lower than the SPC/E, TIP4P 
and TIP5P water models;^ though closer to the experi- 
mental valuej^i^ in accordance with the fact that TIP3P 
water is less structured than the other water models.— 
One may also predict the transition by computing the 
shear viscosity 77. In equilibrium, one can evaluate ry ei- 
ther from the mean square displacement of the Helfand 
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FIG. 10: Temperature dependence of the (a) isobaric specific 
heat, (b) isothermal compressibility. The squares represent 
the experimental values at room temperature.— 

moment associated to viscosity,^^ or from the Green- 
Kubo relation, 

V f°° 

^ = r^/ dt{p^pit)p^pm, (11) 



where V is the volume of the system, T the temperature, 
ks the Boltzmann constant, and Pq^ the pressure tensor, 

= 7 ( E ™»<^f + E I • (12) 

where a and (3 take the values x,y and z. The off-diagonal 
components of permit to compute the shear viscos- 
ity rj while the diagonal ones lead to the bulk viscosity. 
Note that the value of the viscosity calculated by Eq. [11] 
is limited by the length of the simulations. One can then 
estimate the glass transition temperature Tg by fitting 
the power law r] ^ \T — Tg\~^ . We depict in Fig. [TDthe 
temperature dependence of the shear viscosity and the 
power law fit. The fit yields a glass transition tempera- 
ture Tg = 155 K in accord with the values obtained from 
the susceptibilities. 
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